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ABSTRACT 

Aims. We aim to estimate the rate of compact binaries and compare these with similar estimates based on either the observation 
of these binaries or population synthesis models. Since the merger of these compact objects are likely to produce short gamma-ray 
bursts, conclusions about the opening angle of these bursts can be made. 

Methods. We use a set of observed redshift measurements of short gamma-ray bursts to model the rate of these merger events in the 
nearby universe. Various corrections are included in the calculation, such as the field-of-view of the satellite missions, the beaming 
factors of gamma-ray burst and other parameters. 

Results. The computed rate estimates are compared to other rate estimates, based on observations of binary neutron stars and popu- 
lation synthesis models. Given the upper limit established by LIGO/Virgo measurements, it is possible to draw conclusions about the 
beaming angle of gamma-ray bursts. 
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1. Introduction 

Short gamma-ray bursts (GRB) are powerful explosions in the 
depth of the universe, probably created by the merger of two 
compact objects, such as two neutron stars (NS-NS) or a neu- 
tron star and a black hole (NS-BH). These events will be sources 
of strong gravitational waves (GW), which are being searched 
for with ground-based GW detectors such as LIGO and Virgo. 
Efforts have been made to estimate the local rates of these events 
based on the observations of a few known neutron-star binary 
systems and star population models. If short GRBs are indeed 
created by a merger of two compact objects, the observation 
of the redshift distribution can be used to determine the local 
merger rate in a complementary way, which is the main sub- 
ject of this work. Furthermore, known limits on the rate of these 
events from LIGO/Virgo observations can be used to constrain 
parameters related to GRB physics, such as the opening angle of 
the outflows in GRB. 

This paper is organized as follows. After the introduction, the 
available redshift data for short GRBs is reviewed, before the 
fit model is explained and the sources of uncertainties are dis- 
cussed. The results of the calculations are then discussed, and 
compared with other rate estimates and the implications for GRB 
parameters. 



1.1. Gamma-ray bursts 

Gamma ray bursts are intensive bursts of high-energy gamma 
rays, distributed uniformly over the sky, lasting milliseconds to 
hundreds of seconds. Several thousands bursts have been dis- 
covered to date, with the very prominent feature of a bimodal 
distribution of the burst durations, with a minimum around two 
seconds ( |Kouve liotou et al. 1993; Horvath 2002). Bursts with a 



duration shorter than two seconds are called short GRB's, and 
bursts lasting longer than two seconds are labeled long GRBs. 
Long GRBs have been associated with star-forming galaxies and 
core-collapse supernovae ( Campana et al. 2006 Lee et al.| 2004 



Hjor th et al.| [2003; Fruchter 2006[ |Woosley & Bloom| f2006) 



Short GRBs, on the other hand, had been a mystery for a long 
time, but early observations with SWIFT (Gehr els et al.| [2005; 
Barthelmy et al. 2005a) indicated that the progenitor of these 
events is the merger of two compact objects (Eic hler et al.| 1989; 
Narayanetal.|1992| ). 



Besides the merger scenario, some of the short GRBs are prob- 
ably caused by soft gamma repeaters (SGR), rapidly rotat- 
ing magnetically-powered neutron stars, creating sporadic 'star 
quakes' in the crust that generate bursts of gamma radiation 
( |Mereghetti|[2008l | Woods & Thompson|[2004l ). It has been es- 
timated that up to ~25 % of all SG RBs are caused by SGRs 
dTanvir et al.|2005l|Levan et al.|2008| ). 



Throughout this paper, we will use a new GRB classification 
scheme suggested by Zhang et.cd. (Zh ang et al.|2 009), which la- 
bels a GRB with a probable merger progenitor as type-I GRB 
(i.e. the 'short' -GRB class) and GRBs that might have been 
created by a core collapse supernovae as type-II GRB (i.e. the 
Tong'-GRB class). Although this classification is not unambigu- 
ous (there are unclear cases such as GRB 060614, or a sub-class 
might exist related to e.g. soft gamma-repeaters (Chapman et al. 
2008a Czerny et al.|201 1| ), this notation is based on the proba- 
ble underlying physics instead of a single observational quantity, 
the duration of the burst. 
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Table 1. 28 type-I gamma ray bursts with possible redshift measurements and Tqq duration equal to or shorter than four seconds . 
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Notes. The data for these GRB's are taken from (O'Shaughnessy et al. 2008b; Sakamoto et al. 20071 and from diverse GCN circulars 
(http : / /gen . gsf c . nasa . gov/ gcn3_archive . html) . Reliable redshift values are marked in bold, while implausible redshift values 
are marked in italics. ^ The duration refers to the earth frame. No cosmological correction was applied because of the a priori unknown redshift. 



1.2. Associated gravitational wave observations 

As mentioned just above, type-I GRBs might be sources of grav- 
itational waves searched for with gravitational wave detectors, 
such as LIGO and Virgo. The LIGO detectors, described in de- 
tail in ( |Abbott et al.|2004b||Barish & Weiss|1999| l, consist of two 
instruments at two sites in the US, while Virgo consists of one 
instrument located near Pisa in Italy ( Acerne se et al.|2006) l. 
Several results of searches for merger signals have been pub- 
lished, with upper limits on the rates of merger events ( |Abbott 
et al.|2004a||2005a|b|[2Tj06l[2rj07l ), along with trigg ered searches 
for merger signals associated with type-I GRBs ( |Abbott et al. 
2008) |Abadie et aLl|2TJT0] i, but so far, no GWs have been de- 
tected. These detectors are currently undergoing significant en- 
hancements, before a new data-taking campaign commences in 
^2015 with a 10-fold increase in sensitivity with respect to the 
initial configuration, probing a 1000-fold volume in space. 



- Reliable redshifts: GRBs with a very high probability that 
the associated redshift is correct, which is the case when an 
optical counterpart has been identified, or when only one 
galaxy is located in the error-circle of the observation. These 
GRBs are used in datasets labeled 'A' and 'B'. 

- Probable redshifts: GRBs with a good chance that the red- 
shift association is valid, in the case of e.g. two galaxies re- 
maining in the error-circle, or when the observations lead to 
controversial results. These GRBs are used in the dataset la- 
beled 'B'. 

- Implausible redshifts: GRBs with a very low probability 
that the assigned redshift is correct, as in the cases without 
identification of an optical counterpart or when many galax- 
ies reside in the error circle. These GRBs are not used in this 
analysis. 



2. Used data and redshifts 

Because the redshifts are a very crucial piece of information in 
this work, we re-examine the list of type-I GRBs associated with 
a redshift measurement in detail, to determine whether the asso- 
ciation is justified. Table [T] lists all 28 type-I GRBs with a pos- 
sible associated redshift over the five year period that we con- 



sider; the values are taken from (O'Shaughnessy et al. 2008b 
Sakamoto et al. 2007 1 and diverse GCN circulars^ Each redshift 



association is classified into three groups: 



http ://gcn.gsfc.nasa. gov /gen 3 -archive . html 



The full discussion of the classification is given in Appendix |A| 
and summarized in Table [T] In total, 15 GRBs have been clas- 
sified with a reliable redshift association (dataset 'A') and 7 ad- 
ditional GRBs with a probable redshift association (dataset 'B' 
with 22 data-points in total). The 6 GRBs with an implausible as- 
sociation are not considered in this analysis. A bias might exist 
whereby we tend to favor redshift estimates from nearby GRBs, 
since these are in general brighter and the underlying galaxy can 
be identified more easily, but this uncertainty is not taken into 
account. 
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3. Description of the fit model 



3.2. Model-dependent parameters 



We describe the model used to fit the observed redshift values, 
and the parameters needed to convert the fit results into an astro- 
physical rate. 

3.1. The fit model 

We now describe the model used for the fit, which is a gen- 
eral one of the rate of astronomical objects as a function of 
their cosmological redshifts. These models are described e.g. by 
Chapman efaX] ( |20Q8b|a| > and |Guetta & Prraiil ( p005) and are 
used to model the distribution of type-I GRBs. This model ex- 
presses the number of observable objects with a redshift lower 
than some redshift z*, and is given by: 



Several corrections must be applied to the local rate to obtain the 
true rate of binary mergers. These corrections include the beam- 
ing factor of GRBs fr , the field-of-view of the satellite v, the 
fraction of mergers producing a type-I GRB, 77, and the fraction 
of type-I GRBs created by a merger, a. The general expression 
to obtain the true, corrected rate from the uncorrected, local rate 



is 



rcOlT 



r local " 



rj v 



(4) 



R(z) dV(z) 
1 + z dz 



3.2.1 . Beaming factor 

There is convincing evidence that the outflows of GRBs are 



&(L)dL (1) stron gly beamed (Stan ek et al. 1999; H arrison et al. 1999 



In this equation, N(z* ) is the number of type-I GRB above some 
minimum luminosity with a redshift lower than z*, R{z) is the 
rate-function (in units per volume) at a redshift z, Q(L) is the 
luminosity function, and dV(z)/dz is the volume of a comov- 
ing shell at redshift z. The standard cosmology parameters used 
are Hm=0.27, il\=0.13, and h=0.71. The rate function R(z) de- 
scribes the change in the intrinsic rate of objects as a function 
of redshift z, and all functions used for the fits are described 
in Appendix [B] The luminosity function <fr(L) describes the 
distribution of sources as a function of their luminosity, which 
are given in Appendix [C] The upper integration limit is set to 
Lmax = 10 55 ergs, while the lower limit depends on the thresh- 
old of the satellite and the redshift as 



47T 
l + Z 



^LnW^lim 



(2) 



The detection threshold Pu m is taken for the SWIFT satellite, 
and is roughly Pu m ~ 10 _8 ergs/cm 2 as can be seen in Figure 
12 of |Sakamoto etaLl ( |2008] l. 

The fit of the model is performed by a least squares fitting func- 
tion from the scipy module of pythorj^] which uses a modified 
version of the Levenberg-Marquardt algorithm to minimize a 
function representing the difference between the observed data 
and the fitted model. For each dataset, every possible com- 
bination of rate function and luminosity function is used in 
Eq. ([TJ, yielding 16 different fit functions for each dataset. A 
Kolmogorov-Smirnov (KS) probability is calculated to estimate 
the goodness-of-fit of each of these fit functions, and only fits 
with a KS probability of more than 80 % are kept for further 
analysis. 

The norm of this equation, N , denotes the number density of 
GRBs at zero redshift. To obtain the local rate ?'i oca i, one need 
to take into account the observing period (T = 5 years) and the 
fraction of type-I GRBs used for the fit, compared to the total 
number of type-I GRBs observed during this period. The latter 
correction factors are / =28/15 for dataset 'A' and / =28/22 
for dataset 'B'. The local rate is given as 



r local 



[No 
T 



(3) 



This yields the local, uncorrected rate, which has to be corrected 
for several effects, as described in the next subsection. 



http://www.scipy.org 



Grupe et a l. 2006b ." [Burrows et~al . 2006b), which we consider 
in this work. Since gamma-ray bursts are only visible if the 
Earth is inside the cone of the outflow, the true rate will be 
higher than deduced from the redshift fits alone, by a factor 
fr = 1/(1 — cos(#)), the beaming factor. The angle 9 is the 
opening angle of the outflow. 



In the standard model of GRBs (see e.g. Piran (2005 ); Meszaros 
( |2006| l for reviews), the outflowing matter, initially with a 
Lorentz factor of To, is confined to a cone of opening angle 
9. This angle changes when the Lorentz factor of the outflow 
decelerates and becomes comparable to 9 ; the jet then starts 
spreading sideways and an achromatic drop in the light curve 
(jet-break) is expected. Calculations have been performed that 
show a dependency of the jet-break tb on the jet openin g angle 
as 9 = 0.057 (t 3 b /s ( |Sari et al.|[l999| |Frail et al]|200l| i, with ( 
defined as 



c = 



1 



-3/8 



\O.2J 



1/8 / E. 



7.1S0 

10 53 ergs 



-1/8 



0.1cm- 



1/8 



■(5) 



Here z is the redshift of the source, rj^ the efficiency of con- 
verting the energy of the outflow into gamma rays, n is the mean 
circumburst density, and E j so the isotropic equivalent gamma- 
ray energy. 

It should be noted that the observer is in general not directly 
located on the jet axis, but at an angle 9 Q ff relative to this axis 
(van Eerten et al. 2010 1, so that #meas = #tme + # ff- A typical 
observer is more likely to be positioned at a large angle off-axis, 
with a mean value of 



"true 3 
7 off/ = / p(9)9d9 = -9 tme 



(6) 



where we have used the normalized probability density function 
p(9) = 3$ 2 /$true' The true average opening angle of a GRB is 
then given by 



^true — #meas — # ff — -^meas , 



(7) 



which is only about half the angle inferred by measurements of 
the jet-break times as explained above. 

Evidence of jet-breaks in type-I GRBs are very rare, and sum- 
marized in Table [2] The following list gives some more details 
about these findings: 
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Table 2. List of opening angles or their lower limits. 



GRB 


angle in degrees 


beaming factor f b 1 with correction 


050709 


>1.5 a 


< 9000 


050724 


8-15°, >25" 


< 30 or ~90-300 


051221A 


4 C , 3.7 d 


~ 1250 or ~ 1470 


061021 


4.5< = 4.8" 


~870 



Panaitescu 


2006 


i; <W 


Grupe et al. 


2006a i; (c) 


Burrows et al. 


^2006a 


,. id) 


i 1 1 > 

Racusin ct al. 


(2009 



- GRB 050709: T he afterglow of this GRB has was in 

|Panaitescu| ( |2006| l where no jet-break was found within 10 
days of the time of the burst. Two different circumburst den- 
sities have been used (n = 10 cm~ 3 and n — 10 _5 cm -3 ), 
with the latter being more realistic in the case of a type- 
I GRB, yielding an opening angle of 9 > 6° ( Pan aitescu] 
2006). 

- G RB 050724: The af terglow for this GRB was also analyzed 
by |Panaitescu| ( [2006| l who claimed to have seen a jet-break 
one day after the burst time and suggested an opening angle 
of 6 = 10 - 15° or a somewhat tighter value of 6 = 8 - 12°, 
assuming energy injected into the afterglow from long-lived 
X-ray flares (Panaitescu 2006). The density was assumed to 
be in the range n = 0.1 — 1000 cm~ 3 . This result is doubted 
by |Grupe et al.| ( (2006b[ > who claimed that the opening angle 
is larger than 25° for n = 0.1 cm -3 . 

- GRB 051221A: The analysis of the lightcurve for this GRB 
identified three breaks, the last one of which was assumed 
to be the jet break, corresponding to an opening angle of 
9 = 4° or 9 = 8°, depending on the ambien t density 
(n = 10 -4 cm -3 and n = 0.1 cm -3 , respectively) (Burrows 



Dataset A 



et al.|2006b] >. Again, the low-density value seems more ap 



propriate for type-I GRBs. A value of 9 — 3.7° is given in 
Rac usin et all ( |2009) . 
- GRB 061021: The jet-break times and the corresponding 



opening angles for this GRB can be found in |Racusin 



|et al.| ( 2009) , with an assumed circumstellar density of n 
1 cm" 4 . In the case of a type-I GRBs, however, the more 
probable place is a much less dense region, so a choice of 
n~ 10~ 4 cm~ 3 seems more appropriate. Using the mean red- 
shift for type-I GRBs in equation (j5l) and a mean isotropic 
equivalent energy of E 7 jso — 10^° er g> a more realistic 
value for £= 1.0573 has been used. 



The findings of the opening angles of type-I GRBs are summa- 
rized in Table |2j taking into account the correction factor of 4/7. 
The beaming factors range between < 30 to up to ^9000, which 
represents either the observational spread in the real value or 
a multimodal distribution from different processes (i.e. SGR, 
NS-NS merger, NS-BH merger). In the remainder of this paper, 
a value of =500 is used, corresponding to an angle of 
6.3°, which seem to be a reasonable ch oice and agrees with the 
ran ge of 13 < < 10 4 as found in |Soderberg et~alj |2006 1 
and |Levinson et al.| ( |2002| i. We note that these estimates are 
based on only four reliable jet break measurements. This might 
imply that the majority of GRBs have jet breaks at a much later 
time, leading to a much larger beaming factor. However, the 
four presented cases indicate that at least some GRBs have a 
beaming factor in the deduced range. 




00 2.00 4.00 6.00 8.00 10.00 12.00 

Rate [Myr ' Mpc :! ] 



Dataset B 




10.00 20.00 30.00 40.00 50.00 

Rate [Myr ' Mpc :! ] 

Fig. 1. Histogram of the fit results across the different rate func- 



tions and luminosity functions (as listed in Table D.l I. The left 
plot shows the result for dataset 'A' (reliable redshift values), 
while the right plot shows the results for dataset 'B' (reliable and 
probable redshift values). The individual rates are always above 



zero with a minimum rate of 0.25 Myr 1 Mpc 3 for dataset 'A' 



and 0.066 Myr 1 Mpc 3 for dataset 'B'. 



3.2.2. Fraction of type-I GRB originating from a merger 

Not every type-I GRB is caused by a merger event, some possi- 
bly being created from a soft-gamma repeater or different pro- 
cess. As up to 25% of all type-I GRBs might be created by a SGR 
( Tanvir et al.|2005 Levan et al.|2008[ ), an assumption that 75% 
have a merger progenitor seems reasonable; therefore a value of 
a = 0.75 will be assumed in the remainder of this paper. Since 
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rate [L u ; Myr 

750 1000 





0.00 



5.00 



10.00 15.00 20.00 25.00 
rate [Mpc :! Myr '] 



30.00 35.00 



Fig. 2. Ranges of the computed rates using dataset 'A' with only 
the reliable redshift values and dataset 'B' with both the reliable 
and probable redshift values. The red dot indicates the median 
value, the blue bar the 20% to 80% quantile range, while the 
yellow bar covers the total rate range. 



this is the basic assumption on which this work is based, we con- 
sider this value to be fixed. 




1.0 2.0 3.0 5.0 10,0 £0.0 30.0 

Opening angle 9 [degrees] 

Fig. 3. Plot of the rate for model 'delay power', which yields 
the minimum rate (see Appendix IB), as a function of the GRB 
opening angle 9 and the parameter 77. The other parameters are 
kept constant as they are more reliable. The color show the rate 
in units of Myr -1 Mpc -3 , and the black lines delineate the areas 
excluded by GW upper limits in the case of NS-NS and NS-BH 
binaries (LIGO Scientific Collaboration and Virgo Collaboration 
|2010b| l; the lower left areas are excluded. 



3.2.3. Fraction of mergers producing a type-l GRB 

Not every merger will lead to a type-I GRB, because one mate- 
rial object is needed to create the observed ultra-relativistic out- 
flow. This implies that one of the two objects must be a neutron 
star, while the other object is either a neutron star or a black 
hole. Even then, it is hard to estimate the fraction of mergers 
that might create a type-I GRB; a recent investigation inferred 
rj = 0.01 — 0.4 in the case of NS-BH systems pSelczynski et al. 



2007 ), with an even larger spread in the case of NS-NS systems 



( Belczynski et al. 2008 1. The recent discovery of a 2 M Q neutron 
star (Demorest et al. 2010b a) makes it far more plausible that 
two coalescing neutron stars generate type-I GRBs (Ozel et al. 
|2010| >, resulting in a much larger value of r\. A value of 77 = 0.5 
is used as a first guess for the latter analysis, but we investigate 
also the outcome of choosing any possible value between and 
1. 



3.2.4. Field-of-view 

Most of the GRBs considered in this work are detected by 
SWIFT, which has a fiel d-of-view (FOV) of 1.4 sr half-coded 
(Barthelmy et al. 2005b). This corresponds to about 10% visi- 
bility of the sky (v = 0.1) at any time for SWIFT. Although 
GRB data from other missions have been usec0 and the actual 
sky coverage is not constant, a value of 10% seems reasonable 
throughout the period considered in this paper. We also assume 
a 100% duty cycle of the satellites over the entire period of five 
years we consider, which is not exactly true, but the error asso- 
ciated with this parameter is much smaller than the other param- 
eters. 



4. Fit results and discussion 

Following the above discussion of the parameters, the true cor- 
rected rate rcorr can be obtained from the uncorrected local rate 
local w ^ tn ^ e ex P ress i° n 



rcorr = 7500 r local 



(Si 

V500 



f — 
V0.75 



0.50 

V 



0.10 

v 



(8) 



in which the prefactor clearly reflects the default choices (and is 
equal to 500 • 0.75/(0.50 • 0.10)). 



Table D. 1 summarizes the final sample of models, including their 



from HETE, http://space.mit.edu/HETE/ 



fitted parameters and the goodness-of-fit values, when using the 
parameters shown in equation ([8]). For dataset 'A', the lowest rate 
(1.1 Myr -1 Mpc -3 ) is obtained from the model with the delay 
rate function and power luminosity function, while the maxi- 
mum rate (11.3 Myr -1 Mpc -3 ) is obtained from the model with 
a constant rate function and lognormal luminosity function. The 
histogram of the rates are shown, for both datasets, in Figure [T 
and in Figure [2] as a bar plot, indicating the 20%, 50% anc 
80% quantiles. For dataset 'B', the rates are between 2.6 and 
31.6 Myr -1 Mpc -3 . As dataset 'B' includes redshift values that 
are more uncertain and on average closer, it is unsurprising that 
this set yields higher rate values. To be conservative, only results 
from dataset 'A' are used hereafter. 



4. 1. Constraints on GRB parameters 

The effect of the uncertainty in the parameters / b -1 and 77 is 
shown in Fig. [3] which shows the minimum rate obtained from 
a model (i.e. model 'delay power' in Table |D.l) ), as a function of 
the opening angle 9 and the fraction of mergers creating a type-I 
GRB, rj. These lower limits vary between 0.01 Myr -1 Mpc -3 
and 10 4 Myr -1 Mpc -3 , while the upper limits are between 



5 



Alexander Dietz: Estimation of compact binary coalescense rates from short gamma-ray burst redshift measurements 



Table 3. Comparison of realistic rate estimates from this work 
with the estimates given in LIGO Scientific Collaboration and 
Virgo Collaboration! pOTOal l. 





Rlow 


Rre. 


Rhigh 


Rmax 


NS-NS ( Myr -1 Mpc -3 ) 
NS-BH (Myr -1 Mpc -3 ) 
this work ( Myr -1 Mpc -3 ) 


0.01 

6 x 10 -4 
0.03 


l 

0.03 
7.8 


10 

l 


50 

10 4 




Rlow 


Rre 


Rhigh 


Rmax 


NS-NS (L^Myr -1 ) 
NS-BH (L^Myr -1 ) 
this work (L^Myr -1 ) 


0.5 

0.03 

1.5 


50 
1.5 

390 


500 
50 


2500 

5 x 10 5 



Notes. In this table, -Rj ow is identified with the pessimistic estimate (the 
minimum rate obtained in this work), i?pj with the plausible estimate 
(the median rate obtained in this work), i?jjjgjj with the plausible opti- 
mistic estimate, and i?max with the upper limit (the upper limit in this 
work as we ll). The conversion fact or between the two units is 0.0198 
Lio/Mpc 3 (Kopparapu et al. 2008 i. 



0.1 Myr -1 Mpc -3 and 10 5 Myr -1 Mpc -3 . The black lines show 
the 90% upper limits to merger rates as a result of LIGO/Virgo 
searches, which are 43.6 Myr -1 Mpc -3 for NS-BH mergers 
and 172 Myr -1 Mpc -3 for NS-NS mergers (|LIGO Scientific 



Collaboration and Virgo Collaboration 2010b i. The excluded ar- 
eas to the lower left of these lines impose constraints on the 
opening angle 9 of the outflow in GRBs. When assuming that 
most type-I GRBs are created by NS-BH mergers, an opening 
angle of smaller than ~ 1° is excluded. This plot also indicates 
that r\ might not be too small, agreeing with the results obtained 
in[ 



Belczynski eUUT] p007} and |Ozel et aT1 ( [20T0| . 



4.2. Comparison with other rate estimates 

We now compare the rates deduced in this paper with other 
rate estimates. Two cases are considered: The coalescence rate 
of two neutron stars and the rate of coalescence of a neutron 
star with a black hole. For a NS-NS merger the rate is de- 
duced from known binary pulsars in our Milky Way, and was 

although 



-3 



they could be as high as 50 Myr 1 Mpc 3 ( 


LIGO Scientific 


Collaboration and Virgo Collaboration 2010a 


i. The rates pre- 



dieted for NS-BH are much more uncertain, and have been es- 
timated using population synthesis models. Realistic rates are 
around 0.03 Myr " 1 Mpc~ 3 , although they cou ld be as high as 
1 Myr -1 Mpc -3 ( O'Shaughnessy et al.|2008a l. A recent inves- 



tigation of star formation found the rate of NS-BH binaries is 
possibly significantly higher by a factor of ~20, when a lower 
metallicity is assumed in the models ( |Belczynski et al.|2010] >. 
The rates estimated in this work range from 1.1 Myr -1 Mpc - 



to 11.3 



Myr~ 



1 Mpc" 
" 3 for 



with a median value at 
^7.8 Myr -1 Mpc -J for the choice of plausible parame- 
ters used at the beginning of this section. When including 
the uncertainties in the beaming factor fT , ranging from 30 
(corresponding to 26°) to 9000 (corresponding to 1.5°), and for 
the fraction of mergers producing a type-I GRB rj, from 0.01 
to 1.0, the minimum rate becomes 0.03 Myr -1 Mpc -3 , while 
the maximum rate becomes 10200 Myr -1 Mpc -3 . T able 
compares these rate estimates with the ones given in |LIGO 



Scientific Collaboration and Virgo Collaboration (2010a). The 
plausible pessimistic rate estimate Ri ow is identified with the 
minimum rate in the present study, and the realistic rate estimate 
R re is identified with the median value estimated here using the 



plausible parameter choices. This value, ^7.8 Myr -1 Mpc -3 , is 
very sim ilar to the high rate estimate Rhigh for NS-NS mergers 
given in ( |LIGO Scientific Collaboration and Virgo Collabo ration! 
2010a} (see Table [3J. It should also be noted that the maximum 
rate estimate of 10200 Myr -1 Mpc -3 is much higher than the 
excluded values from LIGO/Virgo measurements, which are 
43.6 Myr -1 Mpc -3 f or NS-BH mergers and 172 Myr -1 Mpc -3 
for NS-NS mergers (|LIGO Scientific Collaboration and Virgo 



|Coilaboration|2010b| l. For completeness, Table|3|also shows the 
rates in units of L^ Myr 



5. Summary and conclusion 

This work has utilized redshift measurements of type-I GRBs 
(i.e. short GRBs) to obtain the local rate of NS-NS and NS-BH 
mergers, respectively. The list of available redshifts have been 
revisited in detail, to assess the reliability of each redshift value. 
For the 22 type-I GRBs observed between 2004 and 2009, 15 
have been found to have reliable redshifts, and 7 to have probable 
values; the redshifts of the remaining six type-I GRB were found 
to be too uncertain and were excluded from the analysis. 
A cumulative distribution was constructed using the two datasets 
of redshift values, which had been fitted by models using differ- 
ent functions for the rate and luminosity. A KS-test criterion was 
used to select models with reliable fits. To obtain the true local 
rate, the fit results were corrected for several factors, including 
the fraction of mergers producing a type-I GRB, rj, the number 
of type-I GRBs created by a merger, and the beaming factor of 
the GRBs. 

The obtained rates are consistent with the high-rate estimates 



given in [LIGO Scientific Collaboration and Virgo Collaboration 
(|2010a|), with a median rate of 7.8 Myr -i Mpc -,i . When in- 



Myr -i Mpc" 

eluding the uncertainties from the beaming factor and param- 
eter rj, the rate is found to vary between 0.01 Myr -1 Mpc -3 and 
10 5 Myr -1 Mpc -3 . However, results from LIGO/Virgo observa- 
tions, placing an upper limit on the rates of NS-NS and NS-BH 
mergers, can be used to constrain the opening angles of GRBs; 
the investigation indicates that opening angles of < 1° are ex- 
cluded. 

Further work is required to improve the accuracy of the results, 
i.e. by taking into account the fluxes measured for GRBs, or 
constraining the ranges of some parameters used in this work. 
Ultimately, only direct detections of gravitational waves from 
mergers will yield a more precise rate, and if associated with 
type-I GRBs, the physical properties of gamma-ray bursts will 
also be able to be constrained. 
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Marion and Benoit Mours for carefully reading this manuscript. I also like 
to thank the anonymous referee for useful comments, and the support from 
Astronomy and Astrophysics. 

Appendix A: Detailed redshift discussions 

GRB 020531 HETE detected this short GRB with a du- 
ration of 0.2-1 second [GCN 13990 and follow-up ob- 
servations revealed various sources in the IPN error cir- 
cle [GCN1408,GCN1415], including two new asteroids 
[GCN1400]. One of these sources was found to weaken in 
brightness [GCN1426] very close to an extended object, i.e. the 



4 Each citation from the Gamma-ray Coordinate Netwok (GCN) will 
be given in this format. The citation can be found on the website 

http : / / gen . gsf c . nasa . gov/ gcn3_archive . html 
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probable host galaxy [GCN1427], which has a spectroscopic 
redshift of z=1.0 [GCN1428]. A different fading source was 
found just outside the error box [GCN1434], making this a 
probable redshift estimate only. 

GRB 040924 This is a short GRB with a duration of -1.5 s 
[GCN2754] for which an optical afterglow (OA) has been identi- 
fied [GCN2734]. A spectroscopic analysis of the very likely host 
galaxy [GCN2750] revealed a redshift of z=0.859 [GCN2800]. 

GRB 050416 This is a SWIFT/XRT GRB 
[GCN3264,GCN3268] with a duration of 2.4±0.2 sec- 
onds [GCN3273], and because an optical transient was 
observed [GCN3265,GCN3266] the measured redshift of 
z=0.6535±0.0002 [GCN3542] is very likely. 

GRB 050509B GRB 050509B was found only within the XRT 
error circle [GCN3395] with a radius of 8". Further investiga- 
tions showed that there are at least four more sources in the 
XRT error circle, one of them a probable high-redshift galaxy 
[GCN3401]. The chance association of a low redshift galaxy is 
reported to be very small [GCN3418], hence this redshift esti- 
mation is implausible. 

GRB 050709 Although a long duration is reported for 
this GRB of -130 seconds [GCN3653], the lightcurve 
and other spectral features classify it as a type-I GRB 
[GCN3570,GCN3653]. The association of the afterglow with the 
probable host galaxy [GCN3605,GCN3612] makes the redshift 
estimation of z=0.16 reliable. 



GRB 050724 GRB 050724 has a T 90 duration of strictly 3 ± 1 
seconds, but because it could belong to the type-I GRB class 
[GCN3667] it is considered in the sample. Four objects were 
found in the XRT error circle [GCN3672], of which two are 
identified as Galactic stars [GCN3675,GCN3679]. Additional 
observations hale us improve our confidence, that the object 
labelled as "D" is the host galaxy with a redshift of 0.258 
[GCN3690,GCN3700]. 

GRB 050813 This SWIFT GRB with a duration of 0.6±0.1 
seconds [GCN3793] has been found at a position with several 
faint extended objects, probably forming a galaxy cluster at high 
redshift, which makes this cluster the most likely source of that 
GRB [GCN3798]. Measurements of the redshift of the galaxies 
suggest a value of z=0.722 [GCN3801], although a redshift of 
z=0.65 also seems plausible [GCN3808]. The more conservative 
value of 0.722 will be used in this work as a probable value. 

GRB 050906 SWIFT detected this very short GRB (T 90 = 
0.128 ± 0.016s [GCN3935]) without finding the source with 
XRT, leaving only the BAT error circle to search for an after- 
glow [GCN3927,GCN3935]. Since this circle contains the mas- 
sive star-forming galaxy IC 328 at z=0.031, a galaxy cluster at 
z=0.43, and field galaxies of unknown redshift, any value is un- 
likely. 

GRB 051016B For this SWIFT GRB with a duration of 
4±0.1 s [GCN4104], a probable optical afterglow was found 



[GCN4111], in a galaxy with redshift z=0.9364 [GCN4186]; this 
makes the redshift value reliable. 



GRB 051221 A This GRB had a duration of 1.4±0.2 seconds 
[GCN4365], and an optical afterglow was detected [GCN4375]. 
The measured redshift of z=0.5465 [GCN4384] is therefore re- 
liable. 



GRB 060502B GRB 060502B was a very short GRB 
with a duration of 0.09±0.02 seconds [GCN5064], for 
which two sources were found in the XRT error circle 
[GCN5066,GCN5071]. One is assumed to be a star, while the 
other appears to be an extended object, whose reliable redshift is 
measured to be z=0.287 [GCN5238]. 

GRB 060505 This GRB has a nominal T 90 duration of 4 ± 
lseconds [GCN5142], and therefore not clearly assigned to ei- 
ther type-I or type-II. The position of the optical afterglow was 
found to be 4".3 from the galaxy 2dFGRS S173Z112, with a 
redshift of z=0.089 [GCN5123]. The distance in projection of 
this late-type galaxy was found to be 7 kpc [GCN5123]. No su- 
pernova was found to be associated with this GRB [GCN5161], 
suggesting that this might be either a merger-driven GRB or a 
GRB at a much larger distance. The redshift value is classified 
as probable. 

GRB 060801 This short GRB (duration of 0.5±0.1 s 
[GCN5381]) was found in the SWIFT/XRT instrument 
[GCN5378], which revealed four objects in its field 
[GCN5384,GCN6386]. In the revised XRT error circle 
[GCN5389], two objects remained, of which one is extended. 
The redshift of that extended object is z= 1.131 [GCN5470], 
making it a reliable estimation. 

GRB 061006 An optical afterglow was found for this 
GRB [GCN5718] revealing the r eliable redshift to be 
z=0.4377±0.0002 ( |Berger et dTp007) . The formal duration is 
130±10 s [GCN5704], but initial short spikes lasting ~ 0.5 sec- 
onds, on which SWIFT did not trigger [GCN5702,GCN5710], 
places this GRB in the type-I category. 



GRB 061201 This short GRB (duration of 0.8 ± 0.1 seconds 
[GCN5882]) had an optical afterglow [GCN5896], but no galaxy 
was found at its position. Close-by objects include a galaxy at 
redshift 0.111 [GCN5952], as well as the galaxy cluster Abell 
995, for which a mean redshift of z— 0.0835 was determined 
[GCN5995]. The offset of the GRB from the galaxy would be 
34 kpc in the first case, and 800 kpc in the second case (from the 
center of the cluster). The value of z=0.1 1 1 is used in the further 
analysis^] 

G RB 061 21 This GRB has a nominal duration of T 90 = 85 ± 
5 s [GCN5905], but an initial short spike of duration — 60 ms 
places it into the type-I regime. This GRB has been located in the 
XRT error circle containing three galaxies [GCN5922]. Since no 



Using the value of 0.0835 changes the outcomes of the fits and the 
results of this work insignificantly - it is therefore safe to use a redshift 
value of 0.111. 
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optical transient was found, the association with a given galaxy 
as host is doubtful, which makes the redshift implausible. 



GRB 061217 This short GRB, with a T 90 of 0.30 ± 0.05 
seconds [GCN5930], exhibited no optical afterglow, but the 
XRT position was found to be within 11 arcsec of a galaxy 
[GCN5948]. Two objects have been found in a more deeper 
observation, with the brighter object proposed to be the host 
galaxy of this GRB [GCN5949,GCN5953]. A subsequent obser- 
vation of the proposed host galaxy yields a redshift of z=0.827 
[GCN5965], so the source would have a isotropic -equivalent en- 
ergy release of about 8 x 10 49 erg [GCN5965], which is rather 
high for a type-I GRB. This makes the redshift probable only. 



GRB 070209 The localization of this GRB within the BAT er- 
ror circle contained no single source within the BAT error cir- 
cle, but three X-ray sources in its proximity [GCN6095]. None 
of these sources were found to be the afterglow of this GRB 
[GCN61 19], making the measured redshift of the source closest 
to the GRB position of z=0.314 [GCN6101] very unlikely. 



GRB 071227 For this GRB, an optical afterglow was 
found [GCN7157] to coincide with the single source in 
the XRT error circle [GCN7151]. The redshift of z=0.384 
[GCN7152,GCN7154] is therefore reliable. 

GRB 080121 The only redshift reported for this GRB is 
z=0.046 for two galaxies in the BAT error circle [GCN7210]. 
In this field, many other galaxies are present that might belong 
to a group of galaxies [GCN7210]. Since no XRT position could 
be determined [GCN7209], the redshift value is probable. 

GRB 080520 For this short GRB with a duration of 
2.8±0.7 seconds [GCN7761], an optical afterglow was found 
[GCN7753] for which a redshift of z=1.545 was determined 
[GCN7757]. This is a reliable redshift estimation. 

GRB 09041 7A This GRB with a duration of 0.072±0.018 sec- 
onds [GCN9138], was found to be close to a low-redshift galaxy 
[GCN9134] with a redshift of 0.088 [GCN9136]. Since no op- 
tical afterglow has been found for this GRB, the redshift values 
are implausible. 



GRB 070406 A short GRB with a T 90 of 0.7 ± 0.2 seconds 
[GCN6261], for which a bright source was detected in the XRT 
error circle, whose redshift is 0.703 [GCN6262]. Spectral fea- 
tures indicate this to be a quasar, and unrelated to the burst. 
Further investigations showed no hint of a fading afterglow, and 
the large number of faint sources found in the XRT error circle 
make this redshift estimate of z=0. 1 1 [GCN6249] uncertain. 



GRB 070429B This BAT GRB has a duration of 0.5 ± 0.1 
s [GCN6365], and two faint objects were found in the XRT 
error circle [GCN6372]. For the brighter objects, a redshift of 
z=0.9 023±0.0002 was determined [GCN7104] ( |Cenko et al] 
2008), as well as evidence that this object contains a fading 
source [GCN7145]. Therefore, the measured redshift is reliable. 



GRB 07071 4B This long GRB with a standard T 90 time of 
64±5 seconds shows spectral features of a type-I GRB, espe- 
cially the zero spectral lag [GCN6623]. An optical transient was 
found in the XRT [GCN6630], and the hos t's galaxy redshift i s 
found to be z=0.9225±0.0001 [GCN6836]( jCenko et al.|[2 008 ), 
which makes this a reliable redshift. 



GRB 070724 This SWIFT GRB [GCN6654] had a dura- 
tion of 0.4±0.04 seconds [GCN6656], and two sources were 
found in the XRT error circle, none of which showed variations 
[GCN6673]. The redshift of one of the source was found to be 
z=0.457 [GCN6665], making this a probable estimate. 



GRB 07081 0B This is a short GRB with a duration of only 
0.08±0.01 seconds [GCN6753]. Several possible sources has 
been identified in the XRT error circle, among them a nearby 
bright galaxy at z=0.0385 (source SI) and a cluster of galaxies 
at a redshift of z=0.49 [GCN6756] with an X-ray source (source 
S2 in [GCN6754]). In a latter observation, the second source no 
longer was detected, making this the probable position of the 
afterglow with a redshift of 0.49. 



GRB 090510 For this short GRB with a duration of 
0.3±0. 1 seconds [GCN9337] an optical afterglow has been 
found [GCN9338], which gives a reliable redshift of z=0.903 
[GCN9353]. 

Appendix B: Rate models 

The rate functions used to fit the data according to Eq. ([TJ, except 
for the trivial case of the 'constant' rate, are given by: 

1. The 'sfr' rate that follows the star formation rate (model 
'SF2' ) descr ibed in |Porciani & MadaulpOOT) and |Guetta &] 
|Plranl ( r20031 ): 

R(z) ee R SF2 {z) = R s>0 23 ;! P A (3 ; 4 J 9 • (B.l) 
exp (3.4 z) + 22 

2. The 'merger' rate of two compact objects was derived in 



Guetta & Piran (2005 ) by studying six observed double neu- 
tron stars ( Champion et al.|2004 1. This rate follows a time- 



delay distribution (oc 1 /r) 



R(z) = R M ,o f 
Jo 



dr R S F2 (z(t - t)) jr 



(B.2) 



The 'delay' rate similar to the 'merger' rate, but with a con- 
stant time-delay distribution 



R(z) = R 



Dfl 



t(z) 



dr R SF2 (z(t - T )) 



(B.3) 



Appendix C: Luminosity models 

This section describes the luminosity functions used to fit the 
data according to Eq. ([TJ. 

1 . A single power-law distribution, which is often used to de- 
scribe the pdf of luminosity in astrophysics (with two param- 
eters: $o an d a) 



$(£) - $ 



Lr 



(CI) 
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2. A broken power-law distribution, describing e.g. two un- 
derlying populations in the luminosity ( |Guetta & Piran 2005 \ 
(with four parameters: $ > L , a and p) 



for L < L 



o 



*(L)= $ (£) " fori>=L 



(C.2) 
(C3) 



et al. ( 2Q06) > (with three parameters: <I>o, Lo an d oO 
L 



3. The Schechter distribution, used for example in |Andreon 
;e parameters: 

cxp(-L/L ) 



$(L) = $ 



L 



(C4) 



4. The log-normal distribution, describing a standard candle, 
e.g. a population with about the same luminosity (following 
(Chapman et al.|2008a l, with three parameters: <£>o, Lq, and 



= $ ^ exp 



I /-(logL-logLo) 2 



2a 2 



(C5) 



Appendix D: Fit results 



Table D.l. Fit results with a KS-probability of at least 80 %, 
which have been used in this work. 



Data 


Model 


InLo 


a, a 


KS 


Rate 

[Mpc~ 3 Myr -1 ] 


A 


constant power 


— 


1.8 


0.99 


9.22 


A 


constant schechter 


50.5 


1.8 


1.00 


11.24 


A 


constant lognormal 


24.1 


8.1 


1.00 


11.29 


A 


sfr power 


— 


2.0 


0.81 


7.76 


A 


sfr schechter 


49.3 


1.9 


1.00 


9.20 


A 


sfr lognormal 


45.3 


2.5 


1.00 


8.59 


A 


merger power 


— 


2.2 


0.81 


2.94 


A 


merger schechter 


49.4 


2.2 


1.00 


3.99 


A 


merger lognormal 


44.0 


2.7 


1.00 


3.74 


A 


delay power 


— 


2.3 


0.95 


1.07 


A 


delay schechter 


49.8 


2.3 


1.00 


1.49 


A 


delay lognormal 


35.5 


4.6 


1.00 


1.47 


B 


constant power 




1.9 


0.99 


23.91 


B 


constant schechter 


55.1 


1.9 


0.96 


31.28 


B 


constant lognormal 


-510.2 


37.0 


0.96 


31.25 


B 


sfr power 




2.1 


0.94 


19.12 


B 


sfr schechter 


54.4 


2.2 


1.00 


31.60 


B 


sfr lognormal 


-383.5 


28.7 


1.00 


31.51 


B 


merger power 




2.3 


0.94 


7.04 


B 


merger schechter 


55.3 


2.4 


1.00 


13.30 


B 


merger lognormal 


-317.1 


24.2 


1.00 


13.27 


B 


delay power 




2.3 


1.00 


2.63 


B 


delay schechter 


55.7 


2.5 


0.98 


4.37 


B 


delay lognormal 


-316.0 


23.9 


0.98 


4.35 



9 



Alexander Dietz: Estimation of compact binary coalescense rates from short gamma-ray burst redshift measurements 



References 

Abadie, J. et al. 2010, ApJ, 715, 1453 
Abbott, B. et al. 2004a, Phys. Rev. D, 69, 122001 
Abbott, B. et al. 2004b, Nucl. Instram. Methods, A517, 154 
Abbott, B. et al. 2005a, Phys. Rev. D, 72, 082001 
Abbott, B. et al. 2005b, Phys. Rev. D, 72, 082002 
Abbott, B. et al. 2006, Phys. Rev. D, 73, 062001 
Abbott, B. et al. 2007, Phys. Rev. D, 77, 062002 
Abbott, B. et al. 2008, ApJ, 681, 1419 

Acernese, F. et al. 2006, Classical and Quantum Gravity, 23, S635 
Andreon, S., Cuillandre, J.-C, Puddu, E., & Mellier, Y. 2006, MNRAS, 372, 60 
Barish, B. C. & Weiss, R. 1999, Phys. Today, 52 (Oct), 44 
Barthelmy, S. D. et al. 2005a, Nature, 438, 994 
Barthelmy, S. D. et al. 2005b, Space Sci. Rev., 120, 143 
Belczynski, K., Dominik, M., Bulik, T., et al. 2010, ApJ, 715, L138 
Belczynski, K., OShaughnessy, R., Kalogera, V., et al. 2008, ApJ, 680, L129 
Belczynski, K., Taam, R. E., Rantsiou, E., & van der Sluys, M. 2007, arXiv:astro- 
ph/0703131 

Berger, E., Fox, D. B., Price, P. A., et al. 2007, ApJ, 664, 1000 

Burrows, D. N., Grupe, D., Capalbi, M., et al. 2006a, ApJ, 653, 468 

Burrows, D. N. et al. 2006b, ApJ, 653, 468 

Campana, S. et al. 2006, Nature, 442, 1008 

Cenko, S. B. et al. 2008, arXiv:0802.0874 [astro-ph] 

Champion, D. J. et al. 2004, MNRAS, 350, L61 

Chapman, R., Priddey, R. S., & Tanvir, N. R. 2008a, arXiv:08()2.0()08 [astro-ph] 
Chapman, R., Priddey, R. S., & Tanvir, N. R. 2008b, AIP Conf. Proa, 983, 304 
Czerny, B., Janiuk, A., Cline, D. B., & Otwinowski, S. 2011, New Astron., 16, 

33 

Demorest, P., Pennucci, T., Ransom, S., Roberts, M., & Hessels, J. 2010a, 

arXiv: 1010.5788 [astro-ph.HE] 
Demorest, P. B., Pennucci, T, Ransom, S. M., Roberts, M. S. E., & Hessels, J. 

W. T. 2010b, Nature, 467, 1081 
Eichler, D., Livio, M., Piran, T, & Schramm, D. 1989, Nature, 340, 126 
Frail, D. A., Kulkarni, S. R., Sari, R., et al. 2001, ApJ, 562, L55 
Fruchter, A. S. 2006, Nature, 441, 463 
Gehrels, N. et al. 2005, Nature, 437, 851 

Grupe, D., Burrows, D. N, Patel, S. K., et al. 2006a, ApJ, 653, 462 

Grupe, D. et al. 2006b, ApJ, 653, 462 

Guetta, D. & Piran, T. 2005, arXiv:astro-ph/05 1 1239v2 

Harrison, F. A., Bloom, J. S., Frail, D. A., et al. 1999, ApJ, 523, L121 

Hjorth, J. et al. 2003, Nature, 423, 847 

Horvath, I. 2002, A&A, 392, 791 

Kopparapu, R. K. et al. 2008, ApJ, 675, 1459 

Kouveliotou, C, Meegan, C. A., Fishman, G. J., et al. 1993, ApJ, 413, L101 
Lee, W. H., Ramirez-Ruiz, E., & Page, D. 2004, ApJ, 608, L5 
Levan, A. J. et al. 2008, MNRAS, 384, 541 

Levinson, A., Ofek, E. O., Waxman, E., & Gal- Yam, A. 2002, ApJ, 576, 923 
LIGO Scientific Collaboration and Virgo Collaboration. 2010a, Class. Quant. 
Grav., 27, 173001 

LIGO Scientific Collaboration and Virgo Collaboration. 2010b, 

arXiv: 1005.4655 [gr-qc] 
Mereghetti, S. 2008, Astron. Astrophys. Rev., 15, 225, 0804.025 
Meszaros, P. 2006, Rep. Prog. Phys., 69, 2259 
Narayan, R., Paczynski, & Piran, T. 1992, ApJ, 395, L83 

O'Shaughnessy, R., Kim, C, Kalogera, V., & Belczynski, K. 2008a, ApJ, 672, 
470 

O'Shaughnessy, R. W., Kalogera, V, & Belczynski, K. 2008b, ApJ, 675, 566 
Ozel, F, Psaltis, D., Ransom, S., Demorest, P., & Alford, M. 2010, 

arXiv: 1010.5790 [astro-ph.HE] 
Panaitescu, A. 2006, MNRAS, 367, L42 
Piran, T. 2005, Rev. Mod. Phys., 76, 1143 
Porciani, C. & Madau, P. 2001, ApJ, 548, 522 
Racusin, J. L. et al. 2009, ApJ, 698, 43 
Sakamoto, T. et al. 2007, arXiv:0707.4626v3 [astro-ph] 
Sakamoto, T. et al. 2008, AJ Suppl. Series, 175, 179 
Sari, R., Piran, T., & Halpern, J. P. 1999, ApJ, 519, L17 

Soderberg, A. M., Nakar, E., Berger, E., & Kulkarni, S. R. 2006, ApJ, 638, 930 
Stanek, K. Z., Garnavich, P. M., Kaluzny, J., Pych, W., & Thompson, I. 1999, 

ApJ, 522, L39 
Tanvir, N. R. et al. 2005, Nature, 438, 991 

van Eerten, H., Zhang, W., & MacFadyen, A. 2010, arXiv:1006.5125 [astro- 
ph.HE] 

Woods, P. M. & Thompson, C. 2004, in Compact Stellar X-Ray Sources, ed. 

W. G. H. Lewin & M. van der Klis (Cambridge Univ. Press, Cambridge) 
Woosley, S. E. & Bloom, J. S. 2006, ARA&A, 44, 507 
Zhang, B. et al. 2009, ApJ, 703, 1696 



10 



